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A  theory  is  developed  and  applied  to  determine  whether  an  object  submerged  in  an  ocean  waveguide 
and  insonified  only  by  surface-generated  noise  can  be  detected  with  conventional  sensing  arrays.  An 
expression  for  the  total  noise-field  covariance  of  a  stratified  waveguide  with  a  submerged  object 
present  is  derived  using  full-field  wave  theory.  This  is  evaluated  by  numerical  wave-number 
integration  for  a  spherical  object  in  a  shallow  water  waveguide.  The  Cramer-Rao  lower  bound  on 
detection  error  is  computed  for  several  realistic  shallow  water  scenarios  at  both  low  and  high 
frequency.  The  results  indicate  that  cross-range  localization  is  possible  when  the  array  aperture  is 
sufficient  to  resolve  the  object  scale.  This  conclusion  is  supported  by  beamforming  simulations. 

Range  localization  is  possible  at  greater  distances.  However,  this  requires  high  correlation  between 
direct  and  scattered  noise  fields  at  the  sensor,  which  is  difficult  to  replicate  via  matched  field 
processing.  In  addition,  wave  theory  indicates  that  high  resolution  imaging  of  reflected  ambient 
noise  is  generally  most  effective  within  the  deep  shadow  range  of  the  object.  Beyond  the  deep 
shadow  range,  diffractive  interference  from  the  total  forward  field  may  overwhelm  reflections, 
depending  upon  the  incident  noise  directional  spectrum  and  measurement  range.  Overall,  present 
analysis  indicates  that  the  proposed  detection  scheme  presses  the  limits  of  current  technology. 

PACS  numbers:  43.30.Nb,  43.30. Wi,  43.30.Ft 


INTRODUCTION 

We  investigate  the  possibility  of  detecting  a  submerged 
object  insonified  only  by  surface  noise  in  an  ocean  wave¬ 
guide  at  ranges  which  may  be  much  greater  than  an  object 
dimension.  We  do  so  by  means  of  stochastic  noise  modeling, 
wave  and  scattering  theory,  as  well  as  detection  techniques 
which  are  widely  used  in  the  underwater  acoustics  commu¬ 
nity.  Our  motivation  in  this  effort  is  twofold.  First,  we  are 
interested  in  alternative  methods  for  passive  detection  of  un¬ 
derwater  objects.  Second,  we  are  interested  in  the  ability  of 
submerged  vehicles  and  marine  mammals  to  perceive  their 
environment.  For  example,  can  dolphins  acoustically  image 
prey  without  relying  upon  their  ability  to  actively  insonify  a 
target? 

Our  problem  is  expected  to  have  extremely  low  signal- 
to-noise  ratio  (SNR)  beyond  the  immediate  vicinity  of  the 
object.  This  is  because  the  noise  field  incident  on  an  object  is 
of  the  same  order  as  the  direct  noise  field  measured  at  any 
observation  point.  However,  the  scattered  noise  field  is  re¬ 
duced  by  spreading,  scattering  and  absorption.  If  the  coher¬ 
ence  scale  of  the  noise  at  the  sensor  has  been  augmented  by 
the  object's  presence,  localization  by  coherent  processing 
may  be  possible.  In  order  to  separate  scattered  from  direct 
noise,  which  is  assumed  to  be  azimuthally  homogeneous,  the 
coherence  scale  may  have  to  extend  over  an  aperture  suffi¬ 
cient  to  resolve  the  object.  Detection  is  then  highly  depen¬ 
dent  upon  the  SNR,  array  gain  and  sidelobe  level  of  the 


sensing  device.  Since  these  have  no  simple  geometric  rela¬ 
tionship  to  object  dimension,  range  and  array  aperture  in  a 
waveguide,  the  best  way  to  judge  the  practicality  of  the  con¬ 
cept  is  via  modeling  and  simulation. 

To  our  knowledge,  the  first  documented  discussion  of 
using  naturally  occurring  ambient  noise  in  the  ocean  as  an 
analogy  to  “daylight”  appears  in  a  1985  article.1  Here,  and 
in  a  follow-up  paper,2  the  possibility  of  detecting  submarines 
solely  by  their  noise  absorbing  and  scattering  properties  is 
investigated  by  preliminary  approximations.  More  recent 
experimental,3  and  theoretical,4  results  indicate  the  possibil¬ 
ity  of  detecting  a  submerged  object  using  high  frequency 
ambient  noise.  Extremely  short  range  detections,  within  10 
m  of  a  3-m  object,  were  made  in  the  high  frequency  range  of 
5  to  45  kHz  near  a  beach  with  a  large  pier.  The  experiment 
used  a  parabolic  reflector  directed  at  a  rectangular  object. 
Directional  variations  were  not  reported,  but  only  the  differ¬ 
ence  in  ambient  noise  level  for  differing  orientations  of  the 
object,  reflector  and  beach.  In  the  present  paper,  we  find 
support  for  these  high  frequency  experimental  results  and 
explore  the  possibility  of  detection  at  much  longer  ranges 
and  lower  frequencies  with  full-field  modeling  in  a  realistic 
shallow  water  environment. 

Our  theoretical  formulation  is  exact  for  arbitrary  range 
and  frequency  given  the  exact  scattering  function  for  the  spe¬ 
cific  object  and  environment.  In  order  to  make  the  problem 
more  tractable  for  practical  considerations  we  have  adopted  a 
single  simplifying  approximation.  We  use  the  free-space 
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scattering  function  of  the  object.5  This  ignores  contributions 
from  multiple  scattering  and  is  valid  for  objects  not  too  close 
to  surface  and  bottom  boundaries.  This  approximation  does 
not  significantly  affect  the  accuracy  nor  limit  the  range  or 
frequency  domain  of  interest  for  most  ocean  acoustic  appli¬ 
cations. 

We  adopt  a  fully  stochastic  approach  in  our  theoretical 
formulation,  following  the  noise  modelling  of  Ref.  6.  This 
avoids  biases  and  inaccuracies  associated  with  the  alternative 
approach  of  using  statistical  realizations  to  represent  the 
noise  sources.  The  stochastic  approach  is  also  convenient 
since  the  spatial  covariance  of  the  total  noise  field  on  the 
sensing  array  is  the  quantity  directly  required  for  detection. 
Specifically,  an  object  is  placed  in  the  correlated  noise  field 
of  a  three-dimensional  ocean  waveguide.7  The  noise  field  is 
generated  by  a  continuous  sheet  of  stochastic  sources  below 
the  free  surface  which  extend  uniformly  in  range  and  azi¬ 
muth  from  the  object.  The  azimuthal  directionality  of  the 
noise  field  is  therefore  homogeneous.  (Preliminary  results 
indicate  that  azimuthally  inhomogeneous  noise,  in  the  pres¬ 
ence  of  a  coastline  with  surf,  for  example,  may  be  far  more 
advantageous  to  detection.8  However,  we  reserve  more  thor¬ 
ough  analysis  of  this  for  a  future  article.) 

We  accommodate  sources  near  the  object  in  Sec.  I  by 
modifying  a  previously  developed  theory  for  scattering  from 
an  object  in  a  stratified  waveguide.5  The  former  theory  is  in 
terms  of  normal  modes  and  is  therefore  more  appropriate  for 
large  range  separations  between  noise  source  and  scatterer. 
While  contributions  from  nearby  surface  sources  can  be  in¬ 
corporated  by  artificially  extending  the  basement  layer  to 
approximate  continuum  modes,  we  model  them  precisely 
with  an  alternative  approach.  We  develop  a  spectral  represen¬ 
tation  for  the  scattered  field  using  the  methodology  of  Ref.  5. 
In  Sec.  II,  we  then  integrate  the  scattered  and  direct  field 
contributions  from  the  ensemble  of  stochastic  surface 
sources.  This  leads  to  an  expression  for  the  spatial  covari¬ 
ance  of  the  total  noise  field  in  terms  of  multidimensional 
wave  number  integrals  for  an  arbitrary  object.  These  inte¬ 
grals  are  decomposed  into  the  sum  of  the  scattered,  direct 
and  respective  cross  covariance  contributions  to  the  total 
noise  field.  In  Sec.  Ill,  we  introduce  the  free-space  scattering 
function  of  a  spherical  object.  Spherical  symmetry  enables  us 
to  reduce  these  covariance  terms  to  the  sum  and  product  of 
separate  one-dimensional  wave  number  integrals.  These  ex¬ 
pressions  are  then  evaluated  by  numerical  integration  for 
some  illustrative  examples  using  a  typical  shallow  water 
waveguide.  The  spatial  characteristics  and  coherence  struc¬ 
ture  of  the  noise  is  then  investigated. 

In  Sec.  IV,  we  simulate  detection  for  various  ranges  and 
frequencies  using  conventional  arrays  and  signal  processing 
techniques.  Limitations  on  detection  and  imaging  imposed 
by  diffraction  are  discussed  here  and  in  Appendix  C.  We  then 
compute  the  theoretical  lower  bound  on  object  localization 
error  for  the  given  arrays  using  the  Cramer-Rao  approach. 
This  gives  the  upper  bound  on  performance  and  provides  a 
synoptic  guide  to  some  situations  in  which  detection  by  scat¬ 
tered  surface  noise  is  possible. 


atmosphere 
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FIG.  1.  The  geometry  and  coordinates  of  a  submerged  object  in  a  Pekeris 
waveguide  with  surface  generated  noise.  All  coordinate  systems  have  origin 
at  the  object  center  with  positive  z  pointing  downward.  The  waveguide 
depth  is  h,  the  depth  of  the  object  center  is  d.  The  sound  speed,  density  and 
attenuation  of  the  water  column  and  sediment  are  Cj ,  px ,  a]  and  c2,  p2, 
a2,  respectively. 

I.  FULL-FIELD  SCATTERING  FROM  AN  OBJECT  IN  A 
WAVEGUIDE 

In  our  generalized  stochastic  formulation,  we  model 
sources  of  surface  noise  as  extending  continuously  in  range 
and  azimuth  from  a  submerged  object  in  an  ocean  environ¬ 
ment.  However,  the  expressions  for  the  field  scattered  from 
an  object  in  a  waveguide  of  Ref.  5  are  most  appropriate  for 
large  source  ranges  since  they  are  in  terms  of  normal  modes. 
To  remedy  this,  we  derive  similar  expressions  in  this  section 
which  are  valid  for  arbitrary  source  range  by  reformulating 
the  scattered  field  in  terms  of  wave-number  integrals. 

We  use  the  following  notation:  (;t0,y0,z0)  are  the 
source  coordinates;  (xt,yt,zt)  are  coordinates  on  the  surface 
of  the  object;  and  ( x,y,z )  are  the  field  coordinates.  The  ori¬ 
gin  for  all  of  these  coordinate  systems  is  at  the  center  of  the 
object  with  positive  z  pointing  downward.  Cylindrical  (p,0,z) 
and  spherical  (r,6,<f>)  coordinates  are  defined  as  follows 

xo  =  ro  s*n  cos  <£o>  xt  =  rt  sin  cos 

x  =  r  sin  6  cos  , 

yo=ro  sin  60  sin  <f>0,  yt  =  rt  sin  9t  sin  <f>t , 

y  =  r  sin  6  sin  <f> , 

z0  —  r0  cos  0O,  zt=rt  cos  0t,  z  —  r  cos  6, 

pl=xl+yl  P2t=x2t+yl  P2=x2+y2- 

The  geometry  of  an  object  submerged  in  a  waveguide  with 
surface  generated  noise  is  shown  in  Fig.  1.  For  consistency 
with  the  simulations  to  follow,  the  specified  geometry  is  that 
of  a  shallow  water  Pekeris  waveguide,  although  the  form  of 
the  solution  is  later  generalized  to  include  water  column 
stratification.  The  geometry  of  respective  spatial  and  wave- 
number  coordinate  systems  is  shown  in  Fig.  2. 

We  express  the  total  field  <I>(r)  as  the  sum  of  the  incident 
field  <&,( r)  and  the  scattered  field  <b5(r), 


1704  J.  Acoust.  Soc.  Am.,  Vol.  96,  No.  3,  September  1994 


Makris  et  al.\  Detection  of  a  submerged  object  1704 


Spatial  Coordinates 


Wavenumber  Coordinates 


object 


FIG.  2.  The  geometry  of  spatial  and  wave-number  coordinates.  All  coordinate  systems  have  origin  at  the  object  center. 


<D(r)  =  <I>,(r)  +  <I>s(r),  (1) 

and  begin  with  Eq.  (37)  of  Ref.  5  for  the  scattered  field  from 
a  rigid  surface  with  unspecified  shape 

<I>s(r)  =  —  j  J^[<I,;(rf)  +  <l>J(r()]  G(r\r,)dAt,  (2) 

where  G(r|rf)  is  the  waveguide  Green’s  function,  <&,( r,)  is 
the  field  incident  on  the  object,  and  <£>s(r,)  is  the  scattered 
field  on  the  object. 

We  then  decompose  the  incident  field  into  plane  waves, 

<*>,( r,)~  |  d2^gU,;z,,z0)e-l^Pl-Po)y  (3) 

where  g(£,;z,z0)  satisfies 

+[k(z)2-£]g=~w-  S(z-z0),  (4) 


_P2Ti/(Pi^/)-l 
b'  P2r//(Pi*i)  +  l  ' 


(8b) 


The  vertical  wave-number  magnitude  in  the  sediment  is 

#,=  'J[k2(z)-0],  for  k2(z)^g, 

=  i'I[£-k2{z)],  for  k2{z)<£2,  (9) 

where  k(z)  =  (o/c2  + ia2,  and  c2  and  a2  are  the  sound  speed 
and  attenuation  of  the  sediment.  Substituting  Eq.  (6)  back 
into  Eq.  (3)  we  have  the  field  incident  on  the  object  in  terms 
of  up  and  down  going  plane  waves 


<*>,( r,)=T;  |  d2UA(z0)e‘{-^y,^ 

+  B{z0)ei{-trP<-y^-\e‘SrPo.  (10) 


and 

7f=V[*2Cz)-f?L  for  k2(z)^£2, 

=  i^2-k2(z)],  for  k2(z)<£,  (5) 

is  the  vertical  wave-number  magnitude  in  the  water  column, 
and  k(z)  =  a)/c1  +  ia1,  where  co  is  the  angular  frequency,  c1 
and  ax  are  the  sound  speed  and  attenuation  of  the  water 
column.  After  boundary  conditions  appropriate  to  the  wave¬ 
guide  are  satisfied  (see  Fig.  1)  the  following  wave-number- 
depth  Green’s  function  solution  results  for  h-d^z^z^ 

g(€i;zl,z0)=A(z0)eiy‘z'+B(z0)e~iy>z',  (6) 

where 

1  /  ei7i^0  +  ^)_e-iyiz0\ 

)  <7a) 

and 


Similarly,  we  decompose  the  waveguide  Green’s  function 
into  plane  waves 


G(r|r,)  =  2^  J  (11) 


where 

g(^;znz)  =  C(z)elyZt-\-D(z)e~iyZt ,  for  h  —  d^z^zt, 

=A(z)eiyz<+B(z)e~iyz‘ ,  for  zt>z^-d,  (12) 
and 

1  ieiy(z+2d)-RsRbe2iyhe^iyz\ 


C(z)=- 


4  Triy\  1  ~RsRbe 


liyh 


D{z)  = 


1  {eiT,z-RsRbe2ii'he-i*z+2d)^ 


4  my\  1  -RsRbe 


2  iyh 


(13a) 


(13b) 


-RsiRbie2iy‘h  I eiy‘z°-e-iy‘iz°+2d)\ 

)■ 


(7b) 


Both  factors  are  in  terms  of  the  surface  and  bottom  reflection 
coefficients  Rsi  and  Rbi  for  the  incident  field  where 

*,/=-!,  (8a) 


The  incident  field  subscripts  are  no  longer  appropriate  and 
have  been  dropped. 

Now  both  the  incident  field  on  the  object  ^(r,)  and  the 
waveguide  Green’s  function  G(r|r,)  are  expressed  in  terms  of 
up  and  down  going  plane  waves.  Substituting  these  back  into 
Eq.  (2)  we  have 
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®,(r)  =  -  J  J  e^  Po[A(Zoy(-(r^y^ 

+B(z0)ei(  ~  frP(_  r,z,)]  +  ^(r,)  j 

xA  LL  f  rf2£e'f  P[C(z)e,'(-^'+*2') 

0n,  \27t  J 

+D(2)e,(-«p'-^<)]|dAf.  (14) 

To  incorporate  a  plane-wave  scattering  function,  we  express 
the  wave-number  components  in  terms  of  their  associated 
propagation  angles,  which  may  be  complex.  Once  again,  the 
origin  of  these  coordinate  systems  is  at  the  center  of  the 


object.  A  geometrical  interpretation  of  the  following  expres¬ 
sions  has  been  provided  in  Fig.  2, 

£iX=k  sin  a i  cos  /3t,  £*=&  sin  a  cos  fi, 

£iy  =  k  sin  sin  /?, ,  £y  =  k  sin  a  sin  /3, 

Ji  =  k  cos  at,  y=k  cos  a, 

*?=&+&  e=£+er 

For  example,  p+  yz  =  kr  7}{a,fi;0,(f>),  where 
7}(a,/3;6,<f>)  =  cos  a  cos  0+sin  a  sin  0  cos 

(16) 

is  the  cosine  of  the  angle  between  the  propagation  direction 
(a,0)  and  field  coordinate  direction  (0,<£).  Substituting  these 
into  Eq.  (14)  we  have 


<t>s(r)=-  J  J  |2_  |  d2^  + 

x2-  2-  f  d2g  eiS\C(z)e-ikr>v^-a'p'e>'*<)+D{z)e-ikr>1’(a’P-0<'*‘)~\\dAt, 
dnt  y  2  77  J  ) 


(17) 


for  h  —  d^z^zt.  We  can  then  express  the  scattered  field  in 
the  waveguide  in  terms  of  a  plane-wave  scattering  function 
S (0,<f>',Oi,<f>i)  for  the  object 

(18) 

where 

F(z,z0;i,ii)  =  {A(z0)C(z)S{TT-a,/3;Tr-ai,l3i) 

+A  (z0)D(z)S(  a,/3;  tt- a,  ,/?,) 

+  B(z0)C(z)S(i7- a,j8;a,,/3,) 
+5(z0)Z)(z)S(a,/?;al, /?;)},  (19) 

for  h  —  d^z^Zj.  For  zt>z^  —  d,  substitute  4(z)  for  C(z) 
and  B(z)  for  D(z)  in  Eq.  (19)  above. 

The  above  expression  is  the  exact  solution  for  the  scat¬ 
tered  field  if  S(  0,  <f>\  0,- ,  (fij)  is  the  scattering  function  for  the 
specific  object  with  mixed  boundary  conditions  in  the  spe¬ 
cific  waveguide.  The  expression  is  generalized  to  a  stratified 
water  column  by  reformulating  the  local  planewave  ampli¬ 
tudes  A, B, C,D  to  include  the  appropriate  object  depth  de¬ 
pendence.  However,  the  approximation  that  the  index  of  re¬ 
fraction  is  constant  over  the  object  must  be  valid. 

It  is  difficult  to  obtain  an  exact  expression  for 
S ( 0,  <£;  0i ,  < Pi ) .  Therefore,  we  approximate  S(  6 ,  <f>;  0, ,  <£,)  by 
the  free-space  plane-wave  scattering  function9  S(6,(f> ; 0, , 4 \) 
which  is  most  easily  interpreted  in  the  far  field  where 

(eikr\ 

*rjs(0,^;0,,,^)  (20) 


is  the  scattered  field  at  (0,<£)  for  an  incident  plane  wave 
propagating  in  the  direction  (0,  By  introducing  the  free- 
space  scattering  function,  we  use  the  same  reasoning  as  in 
Sec.  Ill  of  Ref.  5.  The  scattered  field  on  the  surface  of  the 
object  is  approximated  by  the  scattered  field  on  the  surface 
of  the  object  in  free  space.  This  is  valid  so  long  as  multiple 
reflections  between  object  and  waveguide  boundaries  are  not 
important,  i.e.  if  the  object  is  not  too  close  to  boundaries.  If 
we  substitute  S(0,0;07,0f)  for  S(0,<£;  0,  ,<£,)  in  Eq.  (19), 
then  the  resulting  expression  is  exact,  neglecting  multiple 
scattering. 

II.  SURFACE  NOISE  FIELD  COVARIANCE  FOR  A 
STRATIFIED  OCEAN  WAVEGUIDE  WITH  A 
SUBMERGED  OBJECT 

The  spatial  covariance  of  the  surface-generated  noise 
field  is  derived  in  this  section  for  a  stratified  ocean  wave¬ 
guide  with  a  submerged  object  present.  This  is  done  by  com¬ 
bining  the  wave  theory  of  the  previous  section  with  a  sto¬ 
chastic  model  for  the  surface  sound  sources.  We  use  the 
approach  of  Ref.  6,  and  represent  surface  noise  sources  by  a 
continuous  sheet  of  random  monopoles  at  z0  below  the 
pressure-release  surface.  To  simplify  our  derivation,  we  ex¬ 
press  the  total  noise  field  as  a  sum  of  two  parts.  The  first  part 
^(r)  is  the  noise  field  measured  in  the  absence  of  the  ob¬ 
ject,  which  we  refer  to  as  the  direct  noise  field.  The  second 
part  ^5(1*)  is  the  surface  noise  scattered  from  the  object  in 
the  waveguide,  which  we  refer  to  as  the  scattered  noise  field. 
The  spatial  covariance  of  the  total  noise  field  in  the  presence 
of  the  object  is  then 

C(r,r')  =  ([^(r)  +  ^s(r)][^(r')  +  ^(r')]*). 

(21) 
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This  can  be  decomposed  into  the  following  terms: 

C(  r,  r' )  =  CNN(  r,  r ' )  +  CNS(  r,  r' )  +  CSN(  r,  r ' ) 

+  C55(r,r'),  (22) 

where  CNN(r, r')  =  (ty N(r)'P %(*' ))  *s  the  covariance  of 
the  direct  noise  field;  C55(r,r')  =  ('Fs(r)'F|  (r'))  the 
covariance  of  the  scattered  noise  field;  and  CNS(r,r') 
=  (^(r)^l(r'))  is  the  covariance  between  the  direct  and 
scattered  noise  field.  Since  the  total  noise  field  covariance 
must  be  Hermitian,  it  follows  that 


Similarly  we  obtain  the  following  expression  for  the  co- 
variance  between  the  scattered  and  direct  noise  fields: 

CSJV(r,r')  =  (277)2(^-) 

x(i)/l  *2Sd2€  KW’M) 

xg*u\z0;£V[fp+fV1,  (29) 


CW5(r,r')  =  C|iV(r',r).  (23) 

We  first  derive  the  covariance  of  the  scattered  noise  field  by 
integrating  contributions  from  the  infinite  sheet  of  surface 
sources. 


Css(r,r')=J  |  d2r0  d2ro(SJr0)Sl(r'0)) 


X(f>s{  r,r0)<t>*(r',r'0). 


(24) 


where  Sw(r0)  is  the  Fourier  transform  of  the  noise  source 
time  dependence  at  angular  frequency  w  over  the  time  inter¬ 
val  T.  Following  Ref.  6,  we  assume  that  the  spatial  coher¬ 
ence  of  the  noise  sources  depends  only  upon  their  separation 
and  that  they  are  uncorrelated 


<5<u(r0)S:(r')>  = 


2 g2  3flrQ-r0|) 
*2  ko“rol 


(25) 


where  q2!T  is  the  source-strength  spectral  density  of  the 
monopole  sheet.  Substituting  Eq.  (25)  and  Eq.  (18)  into  Eq. 
(24),  and  changing  variables  so  that  pf0  =  R  +  p0,  we  can 
eliminate  the  surface  integral  in  R  to  obtain 

[2  q2 

C55(r,r')  =  27r^-p- 


F(2,20;^fe)F*(2',20;|',^)e4fr-r-r'] 

X  J  d2po  <]-wo. 

Recognizing  that 

J  d2po  e'[fe-i;]-Po=(2lT)2(5(^-^)<$(^-^), 


(26) 


(27) 


we  arrive  at  a  six-dimensional  wave-number  integral 

!2q2\!  1 

k 2 


Css(r,r')  =  (2W)3' 

Xd2$d2%  F(z,z0;£&) 
XF*(z',z 


(28) 


which  cannot  be  simplified  further  without  knowledge  of  the 
scattering  function  of  the  object  and  the  acoustic  properties 
of  the  waveguide. 


where  Cyv5(r,r')  is  available  via  Eq.  (23). 

Finally,  the  direct  noise  field  covariance  is  available 
from  Ref.  6, 


C  nn(  r>r^ '  )  ~  (  2  TT  ) 2 1  -p-j 

x  Vi  di  g(^;z,Zo)g*(^z,,z0Uo(^|p'-p|). 

Jo 

(30) 

We  note  that  for  receivers  at  the  same  depth,  the  direct  noise 
field  covariance  is  real,  and  is  only  a  function  of  horizontal 
separation. 


III.  SURFACE  NOISE  FIELD  COVARIANCE  FOR  A 
STRATIFIED  OCEAN  WAVEGUIDE  WITH  A 
SUBMERGED  SPHERICAL  OBJECT 

Using  the  expressions  developed  in  the  previous  section, 
and  the  relevant  free-space  scattering  function,  we  derive  the 
total  noise  field  covariance  for  a  stratified  ocean  waveguide 
with  a  submerged  spherical  object.  We  choose  a  sphere  prin¬ 
cipally  because  spherical  symmetry  enables  the  multidimen¬ 
sional  integrals  for  the  scattered  noise  field  to  be  reduced  to 
a  combination  of  one-dimensional  wave-number  integrals 
which  are  easily  evaluated  numerically.  In  this  section  we 
present  the  general  solution  and  some  illustrative  examples. 
The  derivation  is  outlined  in  Appendix  A. 

A.  General  solution 

We  introduce  the  following  notation,  which  employs 
standard  symbols  for  Bessel  and  Legendre  functions.  One 
dimensional  wave-number  integrals  with  subscripts  A  or  B 
are  of  the  form 


hA{n,l)=  d^(20)A*(20)P„(y)pf(|), 

(31a) 

3«a(», /,«)=  \jidm^)A^z0)P:[^pf[^. 

(31b) 


Integrals  with  subscripts  C  or  D  are  of  the  form 
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(33a) 


hz;n)=  ( 

Jo 


/CD(p,z;n)=  £d£  C(z)P^-J./0(p£) 


Vp',2';n)-JV  dCMCp't 


r 


+  \_€dS  D{z)Pn\  -j\J0{p£),  (32a) 


XP„  y]A(z0)g*(z\z0;n, 


3 cd(P’z 


,z;n,m)=  ( 

Jo 


£di  C(z)Pmn\-r\Jm(pt) 


d(’ JJ('p’) 


r 


+  |  f  rff  d(z)p;[  j«(pf). 


XP:  y|A(z0)g*(z\z0;f ).  (33b) 


(32b) 


And  integrals  with  subscripts  A  or  B  and  g  are  of  the  form  The  covariance  of  the  scattered  noise  field  is  then 
_ I 


C«(r,r’)-(2„)*(^)(-L)  2  2  /(")r(o|[//M(»./)/cD(P.^«)/cD(p'.*';0 

+ 1  BA(nil)I  Dc(P’z',n)IcD(p'  ,z' ;/)  +  IAB(n,l)ICD(p,z',n)IpC(p'  ,z' ;/)  + I  BB(n,l)I  DC(p,z;n)I^c(p'  ,z'  ',1)] 
mi^")  (»-«)!  (/-«)! 

+  2  2.  (7+^)7  cos  OT(^- 0')[3^(n,/,m)3CD(p,z;n,m)3CD(p,,z';/,/n) 

W=1 

+  3BA(n,/,ffi)3DC(p,z;/i,/n)3^D(p',z';/,wi)  +  3/tB(n,/,m)3CD(p,z;n,m)3^c(p',z';/,m) 


+  3BB(n,/,m)DDC(p,z;n,»j)3DC(p',z';/,OT)]  , 


(34) 


and  the  covariance  between  the  scattered  and  direct  noise  field  is 

vwj_ 

k 2  U* 


C5Af(r,r')  =  (2Tr)4|  -p-|  ( — r]  2  /(«)[  [//tg(p',z';n)/CD(p,z;n)+/Bg(p',z';n)/DC(p,z;n)] 


n  =  0 


\ 


^  {n  —  m)\  \ 

+  2  Zu  (n  +  m)\  cos  m(<l>-  </>')[3Ag(p'  ,z' ;n)3CD(p,z;n,m)  +  3Bg(p'  ,z' ;n)3DC(p,z;n,m)]  , 

m  —  1  / 


(35) 


where  C,vs(r,r')  can  be  obtained  by  Eq.  (23). 

For  a  rigid  sphere  we  have 

j'„{ka) 

/(„,  =  «-  l,.(2„  +  l,-^_,  (36, 

and  for  a  pressure-release  or  soft  sphere 

f(n)  =  i(-m2n  +  l)jp^.  (37) 

Equations  (31)-(35)  are  for  h-d^z**zt.  For  zt>z^—d, 
substitute  A(z)  for  C(z)  and  B(z)  for  D(z)  in  Eqs.  (32a)- 
(32b). 

In  Equations  (34)  and  (35),  the  solutions  are  expressed 
as  a  sum  of  two  sets  of  terms  as  a  result  of  the  addition 
theorem  for  spherical  harmonics.  The  first  set  gives  the  azi- 
muthally  omnidirectional  component  and  is  comprised  of 
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sums  over  zero-order  Bessel  functions  J0  and  Legendre  poly¬ 
nomials  Pn  of  degree  n.  The  second  set  of  terms  gives  the 
azimuthally  directional  component  and  is  comprised  of  sum¬ 
mations  of  higher  order  Bessel  functions  J m  and  the  associ¬ 
ated  higher-order  Legendre  polynomials  P™ .  The  A,  B,  C 
and  D  factors  scale  local  plane-wave  components  of  the 
wave  number-depth  Green’s  function  and  include  informa¬ 
tion  about  waveguide  stratification.  [The  solution  to  the  Pe- 
keris  problem  is  obtained  by  substituting  the  respective  ex¬ 
pressions  for  these  given  in  Eqs.  (7)  and  (13).]  These  factors 
are  weighted  by  the  proper  Legendre  polynomial,  to  give  the 
directional  characteristics  of  the  scatterer,  and  are  converted 
to  field  coordinates  by  a  standard  Bessel  transform  of  the 
above  integrals.  The  asymptotic  form  of  the  spherical  Bessel 
and  Hankel  functions,  and  their  derivatives,  may  be  used  in 
Eqs.  (36)  and  (37)  for  the  large  and  small  ka  cases  respec¬ 
tively  when  the  order  of  the  function  allows  this.  Assymp- 
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totically  large  ka  can  lead  to  smoother  peaks  in  the  wave- 
number  spectrum  and  a  more  stable  numerical  integration 
than  for  a  point  source  at  the  same  frequency. 

The  intensities  of  the  scattered  noise  field,  direct  noise 
field  and  total  noise  field  are  denoted  by  7s(r),  IN{ r),  and 
/( r),  respectively.  These  are  obtained  by  frequency  integra¬ 
tion  of  the  respective  covariances  C55(r,r),  CNN{ r,r),  and 
C(r,r).  For  example,  the  bandlimited  scattered  noise  inten¬ 
sity  about  center  frequency  /0  is 

I  1  Wl\  f/o+A//2 

«r)=dd  CssMdf . 

\PCI\T!  J/0-A//2 

We  adopt  the  convention  that  the  source-strength  spectral 
density  is  normalized  to  unity  such  that  10  log(^2/T)=0  dB 
re:  1  ^tPa2/(Hz  m2).  (Notice  that  //Pam  units  are  consistent 
with  a  standard  monopole  source  strength,  but  in  our  case  are 
also  divided  by  areal  units  for  the  continuous  sheet.)  When 
the  covariances  are  integrated  as  above  over  a  uniform  A / 
=  1-Hz  band,  10  log[pc/5(r)]  has  sound  pressure  level  (SPL) 
units  in  dB  re:  1  pPa. 

B.  Illustrative  examples 

We  present  numerical  simulations  illustrating  the  inten¬ 
sity  of  the  scattered,  direct,  and  total  noise  field  as  well  as  the 
azimuthal  coherence  Rss(p,A<f>)  which  is  defined  in  Appen¬ 
dix  B.  These  examples  give  a  quantitative  description  of  the 
spatial  characteristics  of  the  SNR  as  a  function  of  range  and 
depth  from  the  object  and  the  azimuthal  coherence  scale  of 
the  scattered  noise  field.  This  information  is  important  in 
designing  the  sensor  aperture  and  determining  the  array  gain 
necessary  for  detection. 

All  simulations  are  derived  by  numerical  wave-number 
integration  of  Eqs.  (30),  (34),  and  (35),  for  the  shallow-water 
waveguide  shown  in  Fig.  1.  The  surface  noise  sources  are 
submerged  within  a  quarter  wavelength  of  the  free  surface. 
Each  source  produces  the  dipole  radiation  pattern  consistent 
with  low  to  mid-frequency  field  observations.10  (High- 
frequency  field  observations  are  too  sparse  to  be  included.) 
For  example,  at  low  frequency,  the  radiation  pattern  pro¬ 
duced  by  this  source  placement  can  replicate  that  of  a  bubble 
concentration  near  the  free  surface  which  rapidly  decays  with 
depth.11  This  is  also  the  case  at  higher  frequencies  for  shal¬ 
low  grazing-angle  propagation,  typical  in  a  shallow  wave¬ 
guide,  so  long  as  the  bubbles  remain  acoustically  compact.  In 
either  case,  only  the  absolute  level  needs  to  be  scaled.12 
However,  the  absolute  scale  is  unimportant  for  our  detection 
problem  because  the  scattered  noise  intensity  is  proportional 
to  the  direct  noise  intensity. 

The  noise  sources  insonify  a  sphere  which  is  always  in 
the  center  of  the  waveguide  at  z=0.  The  water  column  and 
sediment  parameters  are  fixed  at  <^  =  1500  m/s,  p^lOOO 
kg/m3,  c2=1700  m/s,  p2=1900  kg/m3,  a2=(0.1/\)  m_1.  Low 
and  high-frequency  examples  are  presented  for  soft  and  rigid 
boundary  conditions  on  the  sphere.  For  the  low-frequency 
cases  the  sources  are  at  z0  =  0 . 5  —  d  meters,  where  h  =  100  m 
and  d=5 0  m.  The  sphere’s  radius  is  a  =  10  m.  A  water  col¬ 
umn  attenuation  of  01=7.OX1O-8  m_1  is  used  for  the  /=50 
Hz  case  where  ka^2;  and  a  water  column  attenuation  of 


a^l^XlO-5  m-1  is  used  for  /=300  Hz  where  ka^  12.6. 
For  the  high-frequency  /=  10  kHz  case,  the  sources  are  at 
z0  =  0. 01875  —  d  meters;  the  sphere’s  radius  is  reduced  to 
a  =0.25  m,  ka**  10.5;  and  the  water  column  attenuation  is 
increased  to  a  x  =0.0003  m-1. 

We  first  consider  the  scattered  noise  intensities  shown  in 
Fig.  3  for  pressure  release  boundary  conditions  on  the 
spheres.  (The  corresponding  scattered  intensities  from  rigid 
spheres  are  not  significantly  different  except  for  the  50-Hz 
case  where  the  overall  levels  are  slightly  lower.)  For  the 
50-Hz  example,  modal  structure  in  the  waveguide  is  evident. 
At  300  Hz,  waveguide  interference  is  less  evident.  For  the 
10-kHz  case,  a  diffuse  beampattern  emerges  from  the  object 
with  asymmetric  directional  variations  about  the  horizontal 
related  to  those  of  the  incident  noise  field.  Here,  the  rapid 
falloff  in  intensity  with  increasing  range  from  the  object  is 
due  to  the  dramatic  decrease  in  a/h  and  \/h.  As  a  result, 
waveguide  interference  is  not  evident  over  the  short  range 
shown.  At  higher  frequency,  the  combination  of  higher  volu¬ 
metric  absorption  in  the  water  column  and  sediment,  as  well 
as  increased  bottom  interaction,  will  significantly  reduce  the 
scattered  intensity  at  increasing  range.  This  is  evident  when 
the  transmission  loss  at  1-km  range  is  compared  for  the  50- 
and  300-Hz  cases,  where  the  sphere  size  remains  constant. 
For  comparison,  we  also  present  the  SPL  due  to  monopole 
point  sources  placed  at  the  center  of  the  waveguide 
r0= (0,0,0),  with  intensity  proportional  to  G(r|r0)G*(r|r0),  in 
Fig.  4.  The  intensity  structure  is  only  similar  for  the  50-Hz 
case  where  ka  **2,  and  monopolelike  behavior  is  anticipated 
in  the  scattered  field  intensity. 

We  next  consider  the  azimuthal  coherence  of  the  scat¬ 
tered  noise  fields  Rss(p,A(f>),  shown  in  Fig.  5.  Note  that  for  a 
monopole  source,  the  phase  fronts  will  be  cylindrically  sym¬ 
metric  and  Rss(p,A<f))=l.  However,  this  is  not  the  case  for 
the  scattered  noise  fields,  unless  ka<  1.  The  azimuthal  coher¬ 
ence  has  distinct  spatial  variations  which  change  with  fre¬ 
quency,  object  size  and  boundary  condition.  For  the  ex¬ 
amples  shown  we  find  that  Rss(p,A<f>)  falls  from  near  unity  at 
azimuthal  separations  A <f)**0  to  less  than  half  this  value  at 
A<f>**\/(2a).  This  demonstrates  that  there  is  an  upper  bound 
on  the  azimuthal  coherence  scale  of  the  scattered  noise.  This 
bound  is  apparently  the  ratio  of  the  only  two  scales  available 
in  the  scattering,  the  wavelength  of  the  noise  \ ,  and  the  di¬ 
ameter  of  the  sphere  2a.  For  future  reference,  we  define  the 
azimuthal  correlation  scale  of  the  scattered  noise  as 
A<f>c^\/(2a).  The  constancy  of  A <f>c  over  increasing  range 
indicates  angular  spreading.  The  fact  that  it  is  a  stationary 
function  of  azimuthal  separation  reflects  the  azimuthal  ho¬ 
mogeneity  of  the  noise.  In  contrast,  the  azimuthal  coherence 
of  a  horizontally  propagating  planewave  of  wavelength  X 
varies  with  horizontal  range  p  about  a  point,  and  is  not  sta¬ 
tionary.  For  example,  the  coherence  peak  is  only  about  the 
propagation  direction,  and  it  decays  to  less  than  half  value  at 
an  angular  increment  on  the  order  of  \l\/p,  when  this  angle 
is  small  enough  to  be  approximately  equal  to  its  sine.  A 
physical  interpretation  for  the  azimuthal  coherence  scale,  as 
well  as  its  ramifications  upon  detection  and  imaging,  is  pre¬ 
sented  in  Appendix  C. 

Finally,  we  compare  the  total  noise  field  intensity  with 
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FIG.  3.  Scattered  noise-field  intensity  for  pressure  release  spheres  centered  in  the  waveguide  with  /i  =  100  m  and  d= 50  m  for  the  frequencies  shown.  Ranges 
extend  from  the  sphere’s  center.  Sphere  radii  are  as  shown.  Sound-pressure  level  is  abbreviated  by  SPL. 
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FIG.  4.  Monopole  source-field  intensity  for  the  waveguide  and  frequencies  of  Fig.  3.  The  point  source  is  centered  at  d=50  m  for  each  example. 
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FIG.  5.  Azimuthal  coherence  for  a  spherical  object  centered  in  the  waveguide  where  /i  =  1 00  m  and  d= 50  m  for  the  frequencies  shown.  The  coherence  is  for 
a  horizontal  plane  passing  through  the  center  of  the  object.  Both  soft  and  rigid  boundary  conditions  are  shown.  Ranges  extend  from  the  sphere’s  center.  The 
azimuthal  correlation  scale  of  the  scattered  noise  is  limited  to  within  \j{2a). 


and  without  the  object  present,  as  is  shown  in  Fig.  6.  The 
total  noise  field  intensity  includes  the  cross  terms  between 
the  scattered  and  direct  noise  field,  see  Eq.  (35).  These  terms 
are  important  since  they  can  be  on  the  order  of  the  scattered 
field  intensity,  see  Eq.  (34),  and  significantly  alter  the  struc¬ 
ture  of  the  total  field  near  the  object.  For  X/h  and  a/h  ap¬ 
proaching  unity,  the  single  scatter  approximation  is  not  valid 
at  very  short  ranges.  This  is  evident  in  Fig.  6  where  the  field 
total  does  not  uniformly  vanish  on  the  a  =  10-m  pressure  re¬ 


lease  spheres  at  50  and  300  Hz.  For  X/h  and  a/h  small,  the 
single  scatter  approximation  is  a  good  one  at  very  short 
ranges.  This  is  evident  in  Fig.  7  where  the  total  field  does 
uniformly  vanish  on  the  0=O.25-m  pressure  release  sphere  at 
10  kHz.  (Situations  in  which  the  cross  terms  may  be  signifi¬ 
cantly  diminished  are  discussed  in  Appendix  D.)  For  the 
higher  ka  cases  of  /=300  Hz  and  /=10  kHz,  omnidirec¬ 
tional  detection  by  scattered  noise  seems  to  be  limited  to 
within  the  deep  shadow  range  of  the  object.  This  range  is 
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1  T  1  I  f  f  I  1  |  1  |  1  |  9  1  I  |  I 

f-30OHz  *-] 
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<-2.0 
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FIG.  6.  Total  noise-field  intensity  with  and  without  a  spherical  object  present  in  the  waveguide  for  h  =  100  m,  d=50  m.  Both  50  and  300  Hz  low  frequency 
cases  are  shown  for  soft  and  rigid  boundary  conditions  on  the  a  =  10-m  radius  sphere.  The  sphere  is  displayed  by  a  black  semicircle  with  white  perimeter. 
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FIG.  7.  Total  noise-field  intensity  with  and  without  a  spherical  object  present  in  the  waveguide  for  high  frequency  /=  10  kHz,  h  =  100  m,  d=50  m.  Both  soft 
and  rigid  boundary  conditions  on  the  a=0.25-m  radius  sphere  are  investigated.  The  sphere  is  displayed  by  a  black  semicircle  with  white  perimeter. 


defined  in  Appendix  C  for  free  space  conditions.  (For  the 
above  frequencies  this  range  translates  to  roughly  18  and 
0.43  m,  respectively.)  For  the  lower  k a  and  higher  \/h  (50 
Hz)  case,  a  shadow  range  is  not  well  defined.  Here,  omnidi¬ 
rectional  detection  is  evidentially  plausible  within  a  water 
depth  in  range. 

IV.  DETECTION 

The  simulations  of  the  previous  section  and  the  analysis 
of  Appendix  C  indicate  that: 

(1)  Object  detection  solely  by  exploitation  of  scattered 
surface  noise  presses  the  limits  of  current  technology.  This  is 
primarily  because  the  SNR  is  typically  near  0  dB  in  the 
vicinity  of  the  object  and  falls  off  rapidly  with  increasing 
range  and  frequency. 

(2)  In  directions  where  the  incident  noise  spectrum 
peaks,  and  ka>  1,  forward  scattering  can  dominate  the  scat¬ 
tered  noise  field.  When  this  occurs,  an  angular  coherence 
scale  of  roughly  \/(2 a)  about  the  object  is  measured,  corre¬ 
sponding  to  the  angular  width  of  the  forward  scatter  peak  for 
a  single  incident  plane  wave.  Under  these  conditions,  in¬ 
creasing  the  aperture  of  a  coherent  array  beyond  this  coher¬ 
ence  scale  will  not  increase  the  array  gain  or  significantly 
enhance  detection  in  a  measurement  of  the  scattered  field; 
and  cross-range  resolution  will  be  limited  to  the  scale  of  the 
object.  The  noise  directional  spectra  for  the  shallow  water 
situations  presented  peak  near  the  horizontal.  This  causes  the 
angular  coherence  scale  just  described  to  occur  in  the  azi¬ 
muthal  plane  about  the  object.  The  azimuthal  homogeneity 
of  the  noise  insures  that  the  coherence  scale  is  stationary  in 
azimuth. 

(3)  High  resolution  imaging  of  reflected  ambient  noise  is 
most  effective  within  the  deep  shadow  range  of  the  object 
a(kal 2)1/3,  defined  in  Appendix  C.  Beyond  this  range  dif¬ 
fractive  interference  from  the  total  forward  field  may  over¬ 
whelm  reflections,  depending  upon  the  incident  noise  direc¬ 
tional  spectrum  and  measurement  range.  (Increasing  object 
size  has  a  proportionally  much  greater  effect  in  extending  the 
shadow  range  than  increasing  frequency.) 

(4)  The  standard  monopole  point  source  field  does  not 
provide  a  good  replica  for  matched  field  processing,  even  for 
a  spherical  object,  unless  ka<  1.  Therefore,  we  avoid 


matched  field  processing  in  this  initial  investigation.  (More 
sophisticated  detection  techniques  are  discussed  in  the  con¬ 
cluding  section.) 

Instead  we  provide  some  examples  of  detection  and  im¬ 
aging  of  the  object  by  plane-wave  beamforming  in  Sec. 
IV  A.  This  is  for  the  array  specified  below  and  the  frequen¬ 
cies  previously  considered.  In  Sec.  IV  B,  we  address  the 
more  general  localization  issue  by  presenting  the  Cramer- 
Rao  bound  on  resolution  of  the  object  in  range,  cross  range, 
and  depth.  We  use  a  soft  sphere  in  all  cases  only  because  the 
corresponding  detection  thresholds  are  slightly  larger. 

For  our  sensor,  we  use  a  simple  upright  billboard  array 
with  N =7X7  hydrophones  at  X/2  spacing  and  length  L  =  (  jN 
—  1 ) A/2.  The  angular  resolution  is  therefore  fixed  for  all 
examples.  The  geometry  of  the  array  with  respect  to  the  ob¬ 
ject  is  shown  in  Fig.  8.  The  beampattern  of  the  array  is 
shown  in  Fig.  9.  A  billboard  array  has  obvious  advantages  in 
its  ability  to  reject  unwanted  noise  in  three  dimensions  as 
well  as  its  ability  to  maintain  relatively  large  array  gain  and 
still  fit  into  a  shallow  waveguide  at  low  frequency. 

A.  Coherent  imaging 

For  detection,  we  use  both  the  standard  Bartlett  proces¬ 
sor 

PBar.(a)  =  At(a)C(a0)A(a),  (38) 


X/2 


FIG.  8.  Detection  geometry  of  the  object  and  sensor  are  shown.  The  sensor 
is  a  7X7  element  upright  billboard  array.  Individual  hydrophones  are  at  X/2 
spacing.  The  broadside  axis  of  the  array  yA  is  in  the  same  vertical  plane  as 
the  sphere’s  center  for  all  examples  in  this  paper. 
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FIG.  9.  The  normalized  beampattern  of  the  billboard  array  described  in  Fig. 
8  steered  to  broadside,  in  the  horizontal  and  vertical.  The  Rayleigh  resolu¬ 
tion  is  X/I  =  19°.  The  main  lobe  occupies  twice  this  angular  width.  The  lobe 
structure  is  that  for  a  rectangular  taper.  Note  that  the  vertical  endfire  direc¬ 
tions  are  projected  onto  a  line  rather  than  a  point  due  to  the  coordinate 
system  chosen  (azimuth  versus  elevation  angle).  The  dB  reference  units  are 
normalized  pressure2/steradian. 


and  the  minimum  variance  distortionless  processor  (MVDP) 
processor,  also  known  as  the  maximum  likelihood  method 
(MLM), 13,14 

^MVDp(a)  =  {At(a)C-1(a0)A(a)}-1.  (39) 

The  normalized  replica  vector  A(a)  contains  the  synthetic 
measurements  A  y  (a)  for  each  sensor  r; .  These  are  evaluated 
at  the  search  coordinate  a.  The  normalized  covariance  matrix 
of  the  total  noise  field 


C(a0)  = 


C(rj,rk) 

2jl,C(r,fr,)’ 


(40) 


is  measured  at  sensors  r and  r*  with  elements  C(iy  ,r*)  from 
Eq  (22).  The  object  is  at  coordinates  a0.  For  example  in 
plane-wave  beamforming,  we  have 


1 

A  ;-(a)  =  —  exp{  —  /  k  •  ( r;  —  rref) } ,  (41) 


where  rref  is  the  location  of  a  reference  sensor  on  the  array. 
Since  the  object  is  sometimes  within  the  Fresnel  zone  of  the 
array,  see  Appendix  C,  we  consider  azimuthally  focused 
plane-wave  beamforming.  In  terms  of  the  array  coordinates 
in  Fig.  8,  the  replica  field  focused  at  broadside  and  range  pA 
is 


Aj(eA,<t>A)=-ex  p 


-ik\  zA.  cos  0A 


A : 


+xA  sin  0A  cos  <bA  —  zr^~ 
Ai  A  VA  2pA 


(42) 


This  includes  a  quadratic  correction  term  in  the  phase.  For 
measurement  ranges  pA<L/2[L/(2\)]1/3,  higher-order  cor¬ 


rections  to  the  phase  are  necessary  otherwise  the  beamformer 
output  will  be  nearly  identical  to  that  of  an  unfocused  array. 
Whereas,  in  matched  field  processing  with  a  point  source,  the 
replica  field  is 

G(r,|a) 

A, (a)  =  - •  (43) 

vS*-iG(r*|a)G*(r*|a) 

We  consider  the  50-Hz  case  first  where  the  vertical  ap¬ 
erture  of  the  array  nearly  spans  the  full  water  column;  ex¬ 
amples  are  shown  in  Fig.  10.  We  begin  with  the  directional¬ 
ity  of  the  noise  as  determined  by  the  MVDP  beamformer. 
Considering  the  vertical  directionality,  two  familiar  peaks  ap¬ 
pear  about  horizontal  broadside  ($4  =  77/2). 14  These  are  asso¬ 
ciated  with  the  discrete  wavenumber  component  of  the  noise 
field.  Steeper  angle  contributions  from  direct  noise  sources 
above  the  array  and  bottom  reflections  are  artificially  in¬ 
creased  by  broadening  at  endfire  in  this  untapered  array.  (The 
coordinate  system  of  the  beamformer  maps  a  single  direction 
to  a  horizontal  line  in  the  rectangular  directionality  plots  for 
both  up  and  down  vertical  endfire.) 

The  directionality  of  the  total  noise  field  in  the  presence 
of  the  spherical  object  at  100-m  range  also  appears  in  Fig. 
10.  The  MVDP  beamformer  is  used  because  it  significantly 
outperforms  performs  the  Bartlett  processor  in  these  ex¬ 
amples.  Although  the  object  is  in  the  array’s  Fresnel  zone, 
including  quadratic  focus  with  respect  to  the  object’s  center 
has  little  effect  on  the  measurement  and  is  not  included  in  the 
example.  A  noticeable  difference  in  noise  directionality  oc¬ 
curs  over  the  azimuth  of  the  object.  This  is  distributed  in 
vertical  angle  due  to  the  waveguide,  i.e.,  the  vertical  angle  of 
the  object  cannot  be  localized  by  plane-wave  beamforming 
due  to  boundary  interference.  A  detection  level  of  roughly 
1.4  dB  is  found  for  an  object  center  at  100-m  range,  and  0.1 
dB  at  0.5-km  range.  Note  from  the  figure  caption  that  detec¬ 
tion  is  significantly  enhanced  by  the  presence  of  cross  terms 
in  the  100-m  detection.  We  consider  the  object  dimension  to 
be  nominally  resolvable  by  the  array  when  the  angle  the 
object  subtends  is  greater  than  or  equal  to  the  \/L  broadside 
(Rayleigh)  resolution.  While  the  object  dimension  is  slightly 
under-resolved  at  100-m  range,  the  example  indicates  that 
detection  is  still  plausible.  At  0.5-km  range,  the  array’s  reso¬ 
lution  is  too  poor  to  resolve  the  object  which  has  become 
practically  undetectable  by  beamforming. 

We  next  consider  the  300-Hz  case,  in  Fig.  11,  for  the 
same  object  dimension,  boundary  conditions  and  ranges 
shown  in  Fig.  10.  This  time  the  frequency-scaled  array  oc¬ 
cupies  15  m  of  the  water  column,  and  the  object  is  in  the  far 
field  of  the  array.  Detection  thresholds  are  similar  to  those  in 
the  50-Hz  examples.  Detection  is  not  significantly  enhanced 
by  the  presence  of  cross  terms.  This  will  be  discussed  further 
in  the  next  subsection. 

Finally,  we  consider  the  10-kHz  case  for  the  reduced 
sphere  diameter  and  much  closer  ranges,  shown  in  Fig.  12. 
The  array  length  is  now  0.45  m,  or  roughly  the  diameter  of 
the  object  which  is  in  the  array’s  far  field.  We  find  a  signifi¬ 
cant  3.6-dB  detection  level  at  2.0-m  range  from  the  object 
center.  For  this  short  range,  bottom  loss  is  insignificant  as  is 
waveguide  interference.  Therefore,  the  object  is  localized  in 
both  range  and  azimuth  by  plane-wave  beamforming.  At  2.0 
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FIG.  10.  Directional  characteristics  of  the  total  noise  field  over  elevation  and  azimuth  angle  are  shown.  These  are  obtained  using  unfocused,  normalized 
MVDP  beamforming  with  the  billboard  array  described  in  Figs.  8-9.  The  dB  reference  units  are  normalized  pressure2/steradian.  The  upper  plots  show  cases 
with  and  without  the  object  present.  These  are  for  /=50  Hz,  /i  =  100  m,  ^=50  m,  and  the  array  center  in  the  middle  of  the  waveguide.  The  object’s  presence 
is  apparent  in  the  beamformer  output.  In  the  lower  examples,  the  beamformer  outputs  with  and  without  the  object  are  differenced  in  decibels  to  highlight  the 
detection  level.  At  100-m  range  from  object  center,  the  detection  level  is  roughly  1.4  dB  and  at  0.5-km  range,  it  is  0.1  dB.  When  covariance  terms  representing 
the  correlation  between  the  direct  and  scattered  noise  fields  are  eliminated,  the  detection  level  is  reduced  to  0.2  dB  for  100-m  range  and  0.05  dB  for  0.5-km 
range. 
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FIG.  11.  Same  as  Fig.  10  except  /=300  Hz.  The  object’s  presence  is  not  apparent  in  the  beamformer  output.  The  object  is  in  the  far  field  of  the  array.  The 
detection  levels  are  similar  to  those  in  the  previous,  lower  frequency,  example.  When  covariance  terms  representing  the  correlation  between  the  direct  and 
scattered  noise  fields  are  eliminated,  the  detection  level  is  1.5  dB  for  100-m  range  and  0.13  dB  for  0.5-km  range.  Their  effect  is  less  significant  than  in  the 
lower  frequency  example. 
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FIG.  12.  Same  as  Fig.  10  except  /=10  kHz,  and  the  object  size  and  boundary  condition  have  been  altered  as  indicated  in  the  figure.  The  object  is  in  the  far 
field  of  the  array.  The  object’s  presence  in  the  beamformed  output  is  apparent  at  2.0-m  range  from  object  center.  The  detection  levels  show  a  maximum  of 
roughly  3.6  dB  at  2.0  m,  and  0.35  dB  at  10.0-m  range.  When  covariance  terms  representing  the  correlation  between  the  direct  and  scattered  noise  fields  are 
eliminated,  these  levels  are  unchanged. 
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FIG.  13.  Cramer-Rao  lower  bound  on  detection  error  for  a  sphere  in  a  waveguide  with  /=50  Hz,  /r  =  100  m,  d=50  m.  The  bound  is  computed  with  and 
without  covariance  cross  terms  included  as  shown.  The  sphere  center  is  in  the  same  horizontal  plane  as  array  broadside.  Only  object-center  ranges  of  50-1000 
m  and  depths  of  -40  to  40  m  with  respect  to  the  array  center  are  shown.  The  array  is  centered  in  the  middle  of  the  waveguide  at  z=0.  The  minimum  detection 
error  for  object  range,  cross  range  and  depth  localization  are  shown. 
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FIG.  14.  Same  as  Fig.  13  except  /=300  Hz. 


m,  the  object  dimension  is  resolvable  by  the  array.  However, 
at  10.0-m  range,  the  object  dimension  is  no  longer  resolvable 
and  detection  is  questionable.  The  effect  of  cross  terms  is 
again  far  less  significant  at  higher  frequency.  In  general,  this 
example  is  not  inconsistent  with  the  high  frequency-short 
range  experimental  results  of  Ref.  3  discussed  in  the  Intro¬ 
duction.  Further,  both  the  array  aperture  and  broadside  reso¬ 


lution  are  roughly  the  same  as  the  respective  aperture  and 
resolution  of  a  dolphin’s  receiving  sonar  at  10  kHz.15  Evi¬ 
dentially,  small  objects  such  as  prey  may  be  passively  detect¬ 
able  by  dolphins  simply  from  their  scattered  noise  fields  at 
high  frequency. 

It  is  noteworthy  that  beyond  the  deep  shadow  range,  the 
10-kHz  scattered-field  intensity  of  Fig.  3  can  now  be  easily 
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FIG.  15.  Same  as  Fig.  13  except /=  10  kHz;  sphere  dimension  and  boundary  conditions  are  as  indicated  in  the  figure;  and  only  object-array  ranges  of  2.0-30.0 
m  and  depths  of  -20  to  20  m  are  shown. 
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interpreted  in  terms  of  forward  scatter  of  the  incident  noise 
directional  spectrum  shown  in  Fig.  12. 


B.  Cramer-Rao  bounds 


We  use  the  Cramer-Rao  upper  bound  on  performance  to 
determine  conditions  in  which  an  object  may  be  detected  by 
scattered  surface  noise.  The  bound  gives  the  optimal  resolu¬ 
tion  obtainable  for  the  given  sensor,  environmental  condi¬ 
tions  and  parameters  to  be  estimated  according  to  estimation 
and  detection  theory.16  We  compute  the  bound  on  range, 
cross-range  and  depth  resolution  of  the  object  for  the  shallow 
water  scenarios  presented  in  previous  examples.  We  begin 
with  Eq.  (32)  of  Ref.  14  which  describes  the  lower  bound  on 
the  mean-square  estimation  error  for  any  unbiased  estimator 
as 


(44a) 


where  elements  of  the  Fisher  information  matrix  are  given  by 


C-’(a) 


<?C(a) 
da ,■ 


C_1(a) 


<?C(a)\ 
9a j  )9 


(44b) 


and  a  is  the  estimated  position,  while  a  is  the  actual  position. 
This  expression  is  valid  under  the  assumption  that  harmonic 
components  of  the  noise  can  be  characterized  by  a  second 
moment,  or  Gaussian  probability  distribution. 

In  all  of  the  examples  presented,  in  Figs.  13-15,  the 
same  X/2  spaced  N =7X7  element  upright  billboard  array  is 
again  placed  in  the  center  of  the  waveguide,  and  oriented 
with  its  broadside  axis  horizontal  and  in  the  vertical  plane 
passing  through  the  object  center.  The  object  is  then  moved 
within  this  plane  in  depth  and  range  from  the  array.  A  plot  of 
the  Cramer-Rao  lower  bound  for  these  respective  object  lo¬ 
cations  is  then  given  on  object  range,  cross  range,  and  depth 
estimation  error.  This  is  done  by  numerically  computing  the 
square  root  of  the  mean  square  estimation  error  given  in  Eq. 
(44a),  where  the  covariances  are  again  determined  by  nu¬ 
merical  wave  number  integration.  To  show  how  correlation 
between  the  direct  and  scattered  noise  field  enhances  detec¬ 
tion,  a  separate  plot  is  given  for  the  root  mean-square  esti¬ 
mation  error  with  and  without  these  covariance  cross  terms 
included. 

Considering  the  50-Hz  case  first,  a  dramatic  reduction  in 
minimum  detection  error  is  found  in  range  and  cross  range 
when  cross  terms  can  be  included.  Also,  when  these  bounds 
are  related  to  the  previous  beamforming  examples,  it  is  clear 
that  the  MVDP  beamformer’s  performance  in  cross  range 
and  depth  resolution  is  near  optimal.  Therefore,  seeking  a 
better  detection  method  for  these  parameters  would  not  be 
necessary.  Also,  it  is  clear  that  the  object  can  only  be  effec¬ 
tively  localized  in  cross  range  and  depth  when  the  angle  it 
subtends  is  resolvable  by  the  array.  This  is  observed  in  the 
previous  beamforming  examples.  On  the  other  hand,  high 
resolution  localization  in  range  is  plausible  out  to  1.0  km.  A 
modified  MFP  technique  may  meet  this  bound,  as  discussed 
in  the  concluding  section. 

At  300  Hz,  the  bounds  are  similar  to  those  found  in  the 
50-Hz  case.  For  example,  cross-range  localization  is  most 
effective  when  the  object  dimension  is  resolvable  by  the  ar¬ 


ray.  However,  beyond  these  relatively  short  ranges,  the  upper 
bound  on  cross-range  resolution  decreases  more  rapidly  for 
the  higher  frequency  case.  This  is  apparently  due  to  the  more 
rapid  decrease  in  SNR  at  higher  frequency  noted  in  previous 
sections.  For  both  50-  and  300-Hz  cases,  when  cross  terms 
are  included,  localization  in  range  consistently  has  the  least 
minimum  error.  This  is  because  the  coherent  structure  of  the 
total  field  is  most  sensitive  to  object  perturbations  in  range. 

Finally,  we  consider  the  10-kHz  case,  which  uses  a 
sphere  of  significantly  reduced  radius.  The  bounds  on  cross 
range  and  depth  again  indicate  that  localization  is  only  plau¬ 
sible  when  the  object  dimension  is  resolvable  by  the  array. 
Cross  terms  only  enhance  detection  for  these  two  parameters 
when  the  array  is  almost  directly  above  and  below  the  object. 
(In  this  case,  the  detection  is  near  endfire  where  front-back 
ambiguity  on  the  array  is  minimized,  but  so  is  resolution.) 
Apparently,  for  these  short-range  orientations,  perturbations 
in  object  cross  range  and  depth  have  as  significant  an  effect 
on  the  field  as  perturbations  in  range. 


V.  SUMMARY  AND  DISCUSSION 

We  have  investigated  the  possibility  of  detecting  a  sub¬ 
merged  object  by  exploiting  the  way  it  disturbs  an  ocean 
waveguide’s  naturally  occurring  surface-generated  noise 
field.  Our  approach  in  this  paper  is  to  formulate  the  problem 
in  terms  of  wave  and  scattering  theory  using  a  stochastic 
model  for  the  noise  sources.  Expressions  for  the  total  noise 
field  covariance  are  derived  for  an  arbitrary  object  in  a  strati¬ 
fied  ocean  waveguide.  Specific  expressions  are  derived  for  a 
spherical  object  and  evaluated  numerically  for  a  shallow  wa¬ 
ter  Pekeris  waveguide.  The  only  simplifying  assumption  is 
that  multiple  scattering  between  the  object  and  the  wave¬ 
guide  boundaries  can  be  neglected.  Detections  are  simulated 
using  a  conventional  upright  billboard  array  design,  with 
7X7  elements  at  X/2  spacing. 

We  have  used  the  Cramer-Rao  upper  bound  on  resolu¬ 
tion  to  determine  some  conditions  in  which  an  object  may  be 
detected  by  scattered  surface  noise  for  several  realistic  shal¬ 
low  water  scenarios.  The  results  indicate  that  (1)  cross-range 
localization  is  possible  when  the  array  aperture  is  sufficient 
to  resolve  the  object  scale.  This  is  supported  by  minimum- 
variance  beamforming  simulations.  Over  ranges  great 
enough  for  bottom  interaction  to  significantly  affect  trans¬ 
mission  of  scattered  noise,  lower  frequency  detections  may 
be  favored  due  to  higher  SNR.  (2)  Range  localization  is  pos¬ 
sible  at  greater  distances.  However,  this  requires  high  corre¬ 
lation  between  direct  and  scattered  noise  fields  at  the  sensor, 
which  is  difficult  to  replicate  via  matched  field  processing 
which  typically  neglects  correlation  between  signal  and 
noise.  (3)  A  simple  monopole  source  replica  is  a  poor  match 
for  the  scattered  field  unless  ka<  1.  A  forward  scattered  wave 
may  prove  a  more  useful  replica.  Since  the  scattered  noise 
field  is  a  stochastic  quantity,  we  believe  that  matched  field 
processing  will  be  of  limited  value  in  meeting  the  Cramer- 
Rao  bound.  It  is  possible  that  a  matched  covariance  proces¬ 
sor  may  prove  to  be  more  optimal,  especially  if  the  replica  is 
obtained  by  full-wave  modeling  of  the  total  noise  covariance, 
as  initial  results  along  these  lines  suggest.17 
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We  have  also  found  that  in  directions  where  the  incident 
noise  spectrum  peaks,  and  ka>  1,  forward  scattering  can 
dominate  the  scattered  noise  field.  When  this  occurs,  an  an¬ 
gular  coherence  scale  of  roughly  \/(2 a)  is  measured  about 
the  object,  corresponding  to  the  angular  width  of  the  forward 
scatter  peak  for  a  single  incident  plane  wave.  Under  these 
conditions,  increasing  the  aperture  of  a  coherent  array  be¬ 
yond  this  coherence  scale  will  not  increase  the  array  gain  or 
significantly  enhance  detection  in  a  measurement  of  the  scat¬ 
tered  field.  These  conditions  also  limit  cross-range  resolution 
to  the  scale  of  the  object  regardless  of  sensor  aperture.  For 
the  shallow  water  scenarios  presented  here,  the  noise  direc¬ 
tional  spectra  peak  near  the  horizontal.  This  causes  the  an¬ 
gular  coherence  scale  just  described  to  occur  in  the  azimuthal 
plane  about  the  object.  The  azimuthal  homogeneity  of  the 
noise  insures  that  the  coherence  scale  is  stationary  in  azi¬ 
muth. 

We  find  that  high  resolution  imaging  of  reflected  ambi¬ 
ent  noise  is  most  effective  within  the  deep  shadow  range  of 
the  object  a(ka!2)xl3 .  Here  the  total  forward  field  will  cause 
the  least  diffractive  interference  in  imaging  reflections.  Be¬ 
yond  this  range,  the  measurement  is  most  likely  to  succeed 
from  an  orientation  such  that  the  incident  noise  directional 
spectrum  does  not  peak  into  the  steering  direction  of  the 
sensor.  Otherwise  the  total  forward  field  may  overwhelm  re¬ 
flections.  However,  maintaining  such  an  orientation  is  diffi¬ 
cult  for  imaging  ranges  greater  than  a  water  depth.  This  is 
because  the  directional  spectrum  of  shallow-water  noise 
tends  to  peak  about  the  horizontal.  For  example,  we  find  that 
including  diffraction  from  the  object  may  substantially  alter 
the  simulated  images  of  Ref.  18.  These  were  obtained  by 
raytrace  assuming  specular  reflection  which  implies  an  infi¬ 
nite  deep  shadow  range.  In  this  reference,  imaging  is  simu¬ 
lated  for  a  sensor  at  20  m  from  a  1-m  radius  sphere  at  75 
kHz.  However,  our  analysis  indicates  that  deep  shadow  ends 
at  5  m.  A  20-m  measurement  is  in  the  range  of  the  Poisson 
diffraction  cone  where  the  total  forward  field  intensity  due  to 
an  incident  plane  wave  is  on  the  order  of  the  incident  field. 
(The  specific  situation  is  discussed  in  more  detail  in  Appen¬ 
dix  C,  where  it  is  shown  that  the  total  forward  field  over¬ 
whelms  reflections  for  incident  plane  waves  of  equal  incident 
intensity.)  Since  75  kHz  is  near  the  upper  frequency  limit  of 
propagating  noise  in  the  ocean,3  this  result  indicates  that 
wave  theory  is  generally  necessary  to  properly  model  acous¬ 
tic  imaging  with  ambient  noise  in  the  ocean.  A  useful  optical 
analogy  might  be  the  theory  for  imaging  extended  non-self- 
luminous  objects  with  incoherent,  coherent  and  partially  co¬ 
herent  light.19  It  is  noteworthy  that  wave  interference  and 
diffraction  caused  by  the  object  is  an  essential  element  of  this 
classic  optical  theory. 

In  summary,  the  detection  scheme  proposed  in  the  title 
of  this  paper  presses  the  limits  of  current  technology  for 
azimuthally  homogeneous  noise.  Recent  work8  indicates  that 
an  azimuthally  directional  noise  field,  such  as  a  coastline 
where  surf  is  generated,  can  provide  a  natural  means  of  sepa¬ 
rating  direct  and  scattered  noise  to  increase  the  detection 
threshold.  In  this  case,  it  may  no  longer  be  necessary  to  have 
aperture  sufficient  to  resolve  the  scale  of  the  object.  Simple 
but  robust  methods  such  as  beamforming  could  be  practical 


at  relatively  long  ranges.  Detections  would  also  be  less  sen¬ 
sitive  to  environmental  fluctuations. 

Finally,  we  note  that  the  array  aperture  and  resolution  of 
our  high  frequency  examples  are  roughly  the  same  as  those 
of  a  dolphin’s  sonar,  operating  in  the  10-kHz  range.15  These 
examples  show  that  a  dolphin  may  be  able  to  detect  small 
objects,  such  as  prey,  within  a  few  meters  range  simply  by 
the  object’s  disruption  of  the  noise  field.  The  Cramer-Rao 
bound  suggests  that  this  is  not  possible  at  much  greater 
ranges  for  the  given  frequency.  The  dolphin’s  evolutionary 
adaptation  of  active  insonification  is  consistent  with  this  re¬ 
sult. 
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APPENDIX  A:  DERIVATION  OF  THE  SCATTERED 
NOISE  FIELD  COVARIANCE  FOR  A  SPHERE  IN  A 
WAVEGUIDE 

We  begin  by  introducing  the  plane-wave  scatter  function 
for  a  sphere 

5(a,yS;a'/,/3,)  =  S(j7(a,yS;a,,)S,)) 


=  2  f(n)P (Al) 
/!  =  () 


Rewriting  the  cosine  of  the  incidence  angle  in  terms  of  the 
vertical  wave  number  we  have 


S( a,/3; a, ,  .ft)  =  5  p-  [  yy-t  +  cos(y8 -  /?,)] 


=  2  /(«) 


«  =  () 


p\z\p  - 

"l  k  "  k 


JL,  (n  —  m)\  y 
+22  t  .  77 p:  - 


",  (n  +  m)\  n\k 


7, 


XP“  -  cosmOS-A) 


(A2) 


where  we  have  used  the  addition  theorem  for  spherical  har¬ 
monics.  We  then  rewrite  Eq.  (28)  as 
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Css(r,r')  =  (277)3(^)(^)  J  j  fd€,d$d?  {l-«'{A(z^0;f,fc)A*(z',z0;f',6)HH 


+  A(z,z0^,^)r*(2',2o;r,^)H  +  -  +  r(z,z0;^^)A*(z',zo;r^,)H-  + 
+  r(z,Zo;|,|I)r*(z',z0;^Mi)H__}, 


(A3) 


where 

A(2,z0;^^)  =  [A(zo)C(z)+B(z0)D(2)], 
r(2,20;^,^)  =  [A(2o)£>(2)+B(2o)C(2)], 


and 


(A4a) 

(A4b) 


2tt 

H±±  =  J  J  |  dpid/3  dp'{e**'-e-P,1 

0 

X5((l/A:2)[±  77,+  ggj  cos (/3— j0,)]) 

XS*((l/*2)[±y'y<+ff,.cos(/8'-A-)])}. 

(A5) 


This  last  expression  can  be  simplified  to  the  following: 


h±±=(2t7)32  2  /(»)/*(/) 

«  =  0  /  =  0 


Mp€)Pn \  ±J, 


XP, 


j\Mp'€')p1 


min  («,/) 

+2  2 

m  =  l 


(w  —  jw)!  (/  —  m)\ 
(n  +  m)\  ( l+m)\ 


cos  m(<f)—  <f>j) 


To  exploit  the  cylindrical  symmetry  of  the  scattered  field 
intensity  and  maintain  a  real  value,  the  azimuthal  coherence 
function  must  be  confined  to  a  horizontal  plane.  For  illustra¬ 
tive  examples  in  this  paper,  the  plane  is  taken  through  the 
middle  of  the  object  and  waveguide  at  zx— 0. 

To  understand  why  the  azimuthal  coherence  can  be  ex¬ 
pressed  by  a  real  function,  it  is  instructive  to  consider  Eqs. 
(34)-(35).  The  covariance  functions  Css( r,r')  and  CSN( r,r'), 
respectively  are  decomposed  into  an  azimuthally  omnidirec¬ 
tional  and  directional  set  of  terms.  This  is  a  consequence  of 
the  addition  theorem  for  spherical  harmonics.  Azimuthal 
variations  in  spatial  coherence  are  determined  by  the 
cos  m (0-0')  factor  in  Eqs.  (34)-(35),  where  m  determines 
the  order  of  the  spherical  harmonics.  Therefore,  for  observa¬ 
tions  points  r  =(p,z,<f>)  and  r'  =  (p,z,<£')  at  the  same  range 
and  depth,  but  differing  azimuth,  the  covariance  of  the  scat¬ 
tered  field  Css(<f>—  <f>')  is  a  real  number  dependent  only  upon 
the  azimuthal  separation,  as  is  CSN((f>—  4>’)  +  CNS(<f>— 
Since  the  intensity  at  a  given  range  is  also  constant,  the  azi¬ 
muthal  coherence  is  a  stationary  function  of  azimuthal  sepa- 
ration  Rss(p;<f>,<f>1)=Rss(p;A<f>)  where  A<£=|<£-<£,|.  This 
is  related  to  the  horizontal  homogeneity  of  the  noise  which 
also  leads  the  direct  noise  field  covariance  to  be  a  real  func¬ 
tion  dependent  only  upon  receiver  separation  for  sensors  at 
the  same  depth. 


xJm(p€)P 


XJm{p'?)P? 


(A6) 


After  substitution  into  Eq.  (A3)  and  some  algebraic  manipu¬ 
lation,  the  various  terms  can  be  separated  into  the  sums  and 
products  of  one-dimensional  wave  number  integrals  shown 
in  Eq.  (34).  After  a  similar  process  the  same  can  be  done  for 
Eq.  (35). 


APPENDIX  B:  AZIMUTHAL  COHERENCE 

The  azimuthal  coherence  can  be  defined  by  a  real,  nor¬ 
malized,  correlation  function 

RSs(P',<t>’<f>i) 

_  Css(P>zi  1&P9Z1  >0i) 

'JCss(p>zi’<f)\  ’P’Zu&^Cssip’Zu&P’Zufi) 

(Bl) 


APPENDIX  C:  ARRAY  DESIGN  AND  IMAGING 
LIMITATIONS  IMPOSED  BY  DIFFRACTION 

We  find  that  diffractive  interference  will  generally  im¬ 
pose  limitations  on  the  ability  to  image  an  object  by  reflected 
ambient  noise  in  the  ocean,  even  at  the  high  frequency  limit 
suggested  in  Ref.  3.  To  demonstrate  this,  we  begin  with  the 
geometry  and  coordinates  of  a  planewave  incident  upon  a 
perfectly  reflecting  sphere  in  Fig.  Cl.  The  incident  wave 
propagates  in  the  <^  =  180°  direction.  Various  regions  about 
the  sphere  are  identified  in  this  adaptation  of  a  diagram  from 
Ref.  9.  Specular  reflection  and  the  limit  of  geometric  rays  is 
a  good  approximation  in  the  lit  region  beyond  the  Fraunhofer 
range  of  the  object,  see  Fig.  C2.  For  shorter  ranges,  higher 
order  corrections  become  appropriate.9  An  exact  expression 
for  the  incident  field  in  the  direction  ^  and  measured  at 
(p,(f>)  is  given  by 

yi  =  eikp  cos(  <£-<£,) 


=  2  in(2n  +  l)j„(kp)P „(cos(<f>—  <f>t)).  (Cl) 

n  =  0 

The  scattered  field  measured  at  (p,<£)  is  given  by 
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Geometrical 
Reflection  Region 


FIG.  Cl.  The  geometry  and  coordinates  of  a  plane  wave  incident  upon  a  perfectly  reflecting  sphere.  The  wave  is  incident  from  the  <p=0  direction.  Various 
regions  about  the  sphere  are  identified  in  this  adaptation  of  a  diagram  from  Ref.  9. 


Vs=-2  »"(2».  +  l)allAj11)(*p)P1I(cos(^-^)),  (C2) 
«  =  0 

where  an=jn(ka)/h[l\ka)  for  a  pressure  release  sphere  of 


Range ,  p 


FIG.  C2.  Graphical  interpretation  of  the  Fraunhofer  condition.  The  angle 
Dfr  subtended  by  an  object  or  array  of  length  D  from  a  distance  p,  where 
p>D .  For  a  point  source  at  p  from  the  center  of  a  broadside  array,  half  the 
aperture  must  be  at  least  p  \jk/p  to  produce  a  path  length  difference  greater 
than  \/4  between  p  and  the  range  from  the  source  to  the  end  of  the  array  for 

p>\. 
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radius  a,  and  the  outward  going  spherical  Hankel  function  is 
used.  This  is  replaced  by  a'n  =  j'„(ka)/h^  ( ka )  for  a  rigid 
sphere. 

The  amplitude  and  phase  of  the  free-space  scattering 
function  S  [Eq.  (Al)]  can  be  derived  from  these 
expressions.9  These  are  plotted  in  Fig.  C3  for  the  ka^  10.5  of 
the  10-kHz  examples  and  a  soft  sphere.  The  peak  of  the 
scattering  function  is  in  the  forward  direction.  This  decays  to 
less  than  half  amplitude  over  an  angular  increment  of 
roughly  \/(2 a).  The  angular  width  of  the  forward  scatter 
peak  roughly  maintains  this  scale  well  into  the  Fresnel  region 
of  the  object.  It  is  no  coincidence  that  this  is  also  the  azi¬ 
muthal  coherence  scale  of  the  scattered  noise  found  in  the 
waveguide  examples  of  Fig.  5.  For  example,  in  directions 
where  the  incident  noise  spectrum  peaks,  and  ka>  1,  forward 
scattering  can  dominate  the  scattered  noise  field.  This  leads 
to  an  angular  coherence  scale  of  roughly  \/(2 a)  about  the 
object,  corresponding  to  the  angular  width  of  the  forward 
scatter  peak  for  a  single  incident  plane  wave.  For  the  shallow 
water  scenarios  presented,  the  noise  directional  spectra  peak 
near  the  horizontal.  This  causes  the  angular  coherence  scale 
just  described  to  occur  in  the  azimuthal  plane  about  the  ob¬ 
ject.  The  azimuthal  homogeneity  of  the  noise  insures  that  the 
coherence  scale  is  stationary  in  azimuth. 

In  coherent  detection,  it  is  desirable  to  have  a  coherent 
signal  over  the  aperture  of  the  sensing  array.  For  example, 
given  N  sensors  at  separations  greater  than  or  equal  to  the 
horizontal  correlation  length  of  the  direct  noise  field  A/2,  the 
coherent  signal  is  elevated  with  respect  to  the  noise  by  the 
array  gain  10  log  N.  Apertures  in  excess  of  the  signal’s  co¬ 
herence  length  scale  are  not  useful  for  detection  since  the 
incoherent  part  of  the  signal  will  be  averaged  out  with  the 
noise.  Therefore,  the  azimuthal  coherence  scale 
A<£c=ss\/(2a)  sets  an  upper  limit  on  the  useful  horizontal 
aperture  for  detection  by  scattered  noise  at  roughly 
Lmax^pA0c,  where  p  is  range  from  sensor  to  object.  Simi¬ 
larly,  the  optimal  azimuthal  resolution  of  the  object  is  limited 
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FIG.  C3.  The  magnitude  (a)  and  phase  (b)  of  the  free-space  plane  wave 
scattering  function  5(0=0°,  <p;  0,=O°,  <ft  =  180°)  for  ka^  10.5  and  a  pressure 
release  sphere.  The  forward  scatter  peak  is  at  <p=180°.  This  peak  decays  to 
half  amplitude  over  an  angular  width  of  roughly  \/(2 a).  The  slope  of  the 
phase  is  on  the  order  of  2 (ka)~l  for  <p  near  90°. 


by  the  Rayleigh  resolution  \/Lm}iX=\/(pA<f>c)  =  2a/ p  to  oc¬ 
cupy  the  angular  scale  subtended  by  the  object.  Under  these 
conditions,  increasing  the  aperture  of  a  coherent  array  be¬ 
yond  this  coherence  scale  will  not  increase  the  array  gain  or 
significantly  enhance  detection  in  a  measurement  of  the  scat¬ 
tered  field.  These  conditions  also  limit  resolution  to  the  scale 
of  the  object  regardless  of  sensor  aperture. 

The  implication  of  these  results  upon  imaging  may  be 
interpreted  in  terms  of  diffraction  by  considering  the  total 
forward  field.  Figure  Cl  also  shows  that  a  diffraction  cone  of 
intensity  comparable  to  that  of  the  incident  field  will  form  in 
the  total  field  beyond  the  deep  shadow.  Therefore,  given  a 
receiver  beyond  the  deep  shadow  range  of  the  object,  the 
total  field  from  a  plane  wave  propagating  in  the  forward 
direction  to  the  receiver  will  dominate  the  field  of  a  reflected 
planewave  because  the  reflected  field  suffers  spreading  loss. 
This  is  confirmed  by  the  results  shown  in  Fig.  C5.  Here  Eqs. 
(C1)-(C2)  are  used  to  compute  the  field  measured  on  a  line 
array  for  the  geometry  shown  in  Fig.  C4.  Figure  C5(a)  shows 
the  magnitude  of  the  field  measured  at  each  sensor.  The  ex¬ 
pected  Poisson  diffraction  pattern  from  the  total  forward  field 
is  evident.  The  peak  of  the  Poisson  diffraction  cone  along  the 
0=180°  axis  is  on  the  order  of  the  incident  field.  The  re¬ 
flected  fields  from  normal  and  back  incidence  plane  waves 
are  orders  of  magnitude  lower.  However,  they  do  not  yet 
follow  inverse  square  spreading.  This  is  found  in  the  Fraun¬ 
hofer  zone  of  the  object  where  the  scattered  field  intensity  in 
the  lit  region  is  equal  to  the  incident  plane-wave  intensity 
divided  by  4 (p/a)2.  Figure  C5(b)  shows  the  measured  direc¬ 
tional  spectra  obtained  by  plane-wave  beamforming  for  the 
same  parameters  as  in  Fig.  C5(a).  The  forward  field  domi¬ 
nates  the  reflected  arrivals  by  about  20  dB  and  shows  peaks 
at  angles  corresponding  to  the  edges  of  the  object  where 
diffraction  begins.  The  central  overlap  corresponds  to  propa¬ 
gation  along  the  Poisson  cone  axis.  The  reflected  arrivals 
have  width  greater  than  the  \/L  =0.9°  broadside  resolution  of 
the  array  because  the  measurement  range  is  well  within  the 
Fraunhofer  range  of  the  array  and  object.  The  direction  of  the 
central  peaks  in  the  two  reflected  arrivals  can  be  associated 
with  specular  reflections  to  the  array  center.  As  expected, 
diffractive  interference  limits  resolution  to  the  scale  of  the 
object.  In  Fig.  C5(c)  all  parameters  are  unchanged  except  the 
array  is  placed  at  p=200  m,  at  roughly  the  Fraunhofer  range 
of  the  array  and  object.  Diffractive  interference  from  the  for¬ 
ward  field  is  greater  than  in  the  previous  example.  The  reso¬ 
lution  of  the  array  is  insufficient  to  distinguish  reflected  ar¬ 
rivals,  which  are  truly  specular.  As  expected,  the  intensity 
ratio  between  the  total  forward  and  specular  arrivals  is  now 
4 (p/a)2.  (Note  the  elevated  sidelobe  levels  well  beyond  the 
azimuth  of  the  object.)  In  Fig.  C5(d)  the  array  aperture,  and 
number  of  elements,  is  increased  by  a  factor  of  3  to  N= 387 
and  replaced  at  p=20  m.  The  main  peaks  of  the  measured 
directional  spectra  show  increased  width.  Increasing  the  ap¬ 
erture  of  the  array  caused  an  increase  in  diverging  beams 
received  from  the  object.  Here  the  array  is  too  close  to  the 
object  for  a  simple  quadratic  focus  to  be  effective,  ie  qua¬ 
dratic  focus  gives  essentially  the  same  result  as  plane-wave 
beamforming  [see  Eq.  (42)].  This  demonstrates  that  an  in¬ 
creased  aperture  without  proper  focus  can  lead  to  poorer 
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FIG.  C4.  Three  incident  plane  waves  of  unit  amplitude  in  /z,Pa  are  scattered  from  a  rigid  sphere  and  received  by  a  line  array  beyond  the  deep  shadow  range 
of  the  object  at  p>ma.  The  array  is  of  length  L  centered  at  <p=180°  with  \/2  spacing.  The  wavelength  of  \=0.02  m  (75  kHz)  and  sphere  radius  a  =  l,0m 
are  taken  from  Ref.  18. 


resolution.  (For  example,  nearsighted  people  can  often  see 
clearly  through  a  pin  hole.)  Since  the  incident  forward  field 
in  the  absence  of  the  object  will  have  measured  intensity  of 
20  log(AT)  for  a  unit  plane  wave,  or  42  dB  for  the  examples 
in  Fig.  C5(b)  and  (c),  arrivals  from  the  object  will  appear  as 
a  silhouette  in  comparison.  (A  simple  optical  example  of  this 
effect  is  to  hold  a  pen  ( a  =5X10-3  m)  at  arms  length  in  front 
of  a  dominant  light  or  window  (\=5X10-7  m).  The  deep 
shadow  ends  at  roughly  ma=0.16  m.  Interference  from  the 
total  forward  field  dominates  specular  reflections  which  are 
not  visible.  When  the  pen  is  brought  to  within  0.16  m  details 
of  the  object  due  to  specular  reflections  from  other  sources 
are  visible.  However,  for  a  much  larger  1-m  radius  object, 
the  deep  shadow  ends  at  a  much  greater  range  of  185  m  in 
daylight!) 

Because  waveguide  noise  directional  spectra  are  typi¬ 
cally  confined  to  two  narrow  bands  about  the  horizontal  (see 
Sec.  IV),  the  total  forward  scattered  field  from  the  object  will 
dominate  reflections  in  the  far  field.  This  is  apparent  in  an 
incoherent  integration  of  S  over  azimuth;  see  Fig.  C3.  Also, 
as  the  resolution  of  the  array  increases,  arrivals  from  the 
object  separate  in  measured  direction  until  the  situation 
shown  in  Fig.  C5(c)  occurs.  In  this  limiting  case,  the  total 
forward  field  from  a  single  incident  plane  wave  is  compared 
to  the  reflected  field  from  a  single  incident  plane  wave.  Here, 
increasingly  large  directional  variations  in  incident  noise,  or 
object  reflectivity,  are  necessary  to  overcome  diffractive  in¬ 
terference  at  increasing  range.  Clearly,  the  measurement  ori¬ 
entation  must  also  favor  reflected  noise  at  the  receiver.  The 
same  conclusion  can  be  drawn  for  an  array  with  a  large  and 
fixed  number  of  receiving  beams  subtending  the  object. 


APPENDIX  D:  DECORRELATION  OF  DIRECT  AND 
SCATTERED  NOISE  FIELDS 

We  consider  two  conditions  in  which  the  covariance  be¬ 
tween  the  direct  and  scattered  noise  field  may  vanish.  The 
first  concerns  relative  motion  between  sources,  scatterer  and 


receiver.  The  second  concerns  temporal  decorrelation  of  the 
noise  source  signal  during  travel  to  the  receiver  from  the 
scattered  and  nonscattered  paths. 

In  the  former,  a  scattering  object  may  under  go  random 
motion  A  r  with  respect  to  the  sound  sources  and/or  receivers 
during  the  sensing  time  window  T.  In  this  case  the  covari¬ 
ance  cross  terms  CSN(ry)+CNS(r9r')  may  become  negligi¬ 
bly  small  when  relative  motion  is  larger  than  a  wavelength. 
For  example,  the  expectation  value  of  the  cross  terms 
(CNS(r9r'  +Ar))Ar  may  vanish.  To  demonstrate  this,  in  Eq. 
(35)  we  replace  ICD(p' ,z' ;n),  IDC(p' ,z' ;n),  3Cb(p'> 
z'\n,m),  and  IDC(p',z';n )  by  (ICD(p'  +  Ar,z' ;n))Ar, 
(hcip’  +  Ar,z';n))  Ar,  0DC(p'  +  Ar,z' ;n,m))Ar,  and 
Ocd(p'  +Ar,z';n,m))frr,  respectively,  where  we  have  as¬ 
sumed  the  change  in  azimuthal  angle  A <t>  is  negligible. 
Then,  for  example, 


OcD(P  +  Ar,z-,n,m))^r 


c(z)p:[^+d{z)p :( 
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FIG.  C5.  Simulated  measurements  for  the  detection  geometry  and  parameters  shown  in  Fig.  C4.  In  all  examples,  the  scattered  field  is  computed  by  the  exact 
summation  of  Eqs.  (C1)-(C2):  (a)  the  magnitude  of  each  of  the  three  fields  measured  on  a  N= 129  element  array  at  p=20  m;  (b)  the  measured  directional 
spectrum  of  the  same  array  as  in  (a).  The  reflected  arrivals  have  width  greater  than  the  \/L  =0.9°  broadside  resolution  of  the  array  because  the  measurement 
is  well  within  the  Fraunhofer  zone  of  the  array  and  object,  (c)  The  AT=129  element  array  is  roughly  at  the  Fraunhofer  range  of  the  array  and  object  at  p= 200 
m.  Forward  scattering  interference  is  greater  than  the  previous  example.  The  resolution  of  the  array  is  insufficient  to  distinguish  reflected  arrivals,  which  are 
truly  specular,  (d)  The  array  aperture,  and  number  of  elements,  is  increased  by  a  factor  of  3  at  p=20  m.  The  main  peaks  of  the  measured  reflected  spectra  show 
increased  width. 
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Here  we  have  used  the  far-field  asymptote  for  the  Bessel 
function,  and  have  assumed  that  p>Ar.  The  expectation  val¬ 
ues  and  therefore  the  cross  terms  C  SN(r,rf)+C  NS(r,rf)  them¬ 
selves  vanish  when  Under  these  circumstances, 

steep  angle  propagation  will  dominate  the  cross  terms.  Fur¬ 
thermore,  due  to  higher-order  mode  stripping,  the  cross  terms 
will  be  most  significant  in  close  proximity  to  the  scattering 
object.  However,  for  the  relative  motion  to  have  any  effect  at 
all,  even  on  perfectly  horizontal  propagation  at  ij=k,  the 
scale  of  the  fluctuations  A  r  must  be  greater  than  X/2.  There¬ 
fore,  this  argument  is  only  of  significance  in  the  mid  to  high 
frequency  regime. 

We  note  that  by  similar  arguments  it  can  be  shown  that 
(C55(r,r'))Ar  does  not  vanish  and  is  relatively  invariant  to 
small  random  object  motions  regardless  of  range  and  propa¬ 
gation  angle.  This  is  due  to  the  presence  of  squared  factors 
such  as  (cos2[Ar(f-£')])A,  and  (sin2[Ar(£-£')])Ar  in  the 
scattered  field  covariance. 

The  covariance  cross  terms  may  also  be  severely  dimin¬ 
ished  by  temporal  decorrelation  of  the  source  signal  over  the 
travel  time  difference  between  source-object-receiver  and 
source-receiver  paths.  In  general,  the  direct  and  scattered 
noise  fields  will  be  independent  if  the  acoustic  travel  time 
between  scatterer  and  receiver  is  greater  than  the  temporal 
correlation  scale  of  the  noise  source.  Since  a  typical  correla¬ 
tion  time  scale  for  a  surface  noise  source  is  on  the  order  of 
r=  1  s,20  we  deduce  that  covariance  cross  terms  will  not  be 
negligible  for  detections  within  roughly  cr/2=750  m  of  the 
object.  Also,  the  contribution  to  covariance  cross  terms  will 
not  be  negligible  from  sources  where  the  path  length  differ¬ 
ence  in  source-object-receiver  versus  source-receiver  paths 
is  less  than  cr/2.  This  condition  selects  sources  in  a  sector 
extending  behind  the  object  from  the  sensor  which  is  sym¬ 
metric  about  the  line  joining  the  two. 
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